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PACS: 24.10.Cn; 24.10.-i; 25.70.-z; 25.75.-q 
Keywords: Many-body theory; Nuclear-reaction models and methods; Low and interme- 
diate energy heavy-ion reactions; Relativistic heavy-ion coUisions 



* supported by GSI Darmstadt 



1 



1 Introduction 



The many-body theory of strongly interacting particles out of equilibrium is a challenging 
problem since a couple of decades. Many approaches based on the Martin-Schwinger hier- 
archy of Green functions have been formulated 0, ^, ^, ^ and applied to model cases. 
Nowadays, the dynamical description of strongly interacting systems out of equilibrium is 
dominantly based on transport theories and efficient numerical recipies have been set up 
for the solution of the coupled channel transport equations ^ ^, 0, [11], |T2|, ^ (and 
Refs. therein). These transport approaches have been derived either from the Kadanoff- 

in Refs. [ITS 



Baym equations in 



O, 18, 19| or from the hierarchy of connected 



equal-time Green functions 0, ^ in Refs. 0, 0] by applying a Wigner transforma- 
tion and restricting to first order in the derivatives of the phase-space variables {X, P) . 

However, as recognized early in these derivations 0, |17|, the on-shell quasiparticle 
limit, that invoked additionally a reduction of the 8iV-dimensional phase-space to 7N 
independent degrees of freedom, where N denotes the number of particles in the sys- 
tem, should not be adequate for particles of short life time and/or high collision rates. 
Therefore, transport formulations for quasiparticles with dynamical spectral functions 
have been presented in the past 0, ^ providing a formal basis for an extension of the 
presently applied transport models §, p|, |T2|, |23|, M, |2|, M, m, m,m,m. Another branch 



of extensions has been developed to include stochastic fiuctuations in the collision terms 
pi] , P^ , P^ ; these models will be denoted as Boltzmann-Langevin (BL) approaches 
in spite of the different numerical approximation schemes. Apart from the transport ex- 



tensions mentioned above there is a further branch discussing the effects from nonlocal 

we will restrict to local 



collisions terms 



^ ^, |3^, p^, p8|, |3^, In this work, however, 

collision terms (for simplicity) and concentrate on the propagation of dynamical spectral 
functions. Recipies to include nonlocal effects in collision terms can be simulated in line 
with Ref. H. 

So far the extension of transport theory to dynamical spectral functions is limited to 
the formal level and only a few attempts have been made to simulate the dynamics of 
broad resonances pT| , In the nuclear physics context this is of particular importance 
for dilepton studies in 'j,n,p + A reactions as well as nucleus-nucleus collisions since 
the vector mesons p and a; are expected to change their properties, i.e. their pole mass 
and width, during the propagation through the nuclear medium P3| , P| , Apart from 
e^e~ in-medium spectroscopy the broadening of the hadron spectral functions in the 
medium should also have some influence on 'subthreshold' meson production; here again 
the question is if an enhanced yield might be due to a downward shift of the meson pole 
mass or simply due to the broadening of its spectral function. Up to now no consistent 
treatment of both, i.e. the real and imaginary parts of the hadron self-energies, is available. 
In this work we particularly address this question and develop a semiclassical transport 
approach for space-time dependent complex self-energies. 

The paper is organized as follows: In Section 2 we will derive the generalized trans- 
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port equations on the basis of the Kadanoff-Baym equations [T^ following in part the 
formulation of Henning [ffS]. Special emphasis will be put on the difference to the con- 
ventional approaches and to the derivation for an energy conserving semi classical limit. 
In Section 3 the resulting set of equations of motion is solved for the case of a complex 
time-dependent Woods-Saxon potential which allows to transparently demonstrate the 
particle propagation in the 8-dimensional phase space of a testparticle. A first applica- 
tion to nucleus-nucleus collisions is presented in Section 4 for energies up to 95 A MeV 
showing in detail the influence of nucleon off-shell propagation on the high energy proton 
spectra and production of high energy 7-rays. The calculational results here are controlled 
by experimental data from GANIL 0. We, furthermore, calculate upper limits for 



kaon production in Ar + Ta reactions at 92 A MeV and compare to the experiment from 
Ref. A summary and discussion of open problems concludes this study in Section 5. 



2 Derivation of semiclassical transport equations for 
particles with dynamical life times 

In this Section we briefly recall the basic equations for Green functions and particle self- 
energies as well as their symmetry properties that will be exploited in the derivation of 
transport equations in the semiclassical limit. 



2.1 Properties of spectral functions 

Within the framework of the closed-time-path formalism 0, ^ Green functions S and 
self-energies S are given as path ordered quantities. They are defined on the time contour 
consisting of two branches from (+) —00 to 00 and (-) from 00 to —00. For convenience 
these propagators and self-energies are transformed into a 2 x 2 matrix representation 
according to their path structure 0, i.e. according to the time branches or (-) for t 
and t'. Explicitly the Green functions are given by 



(1) 



where the subscript -^y denotes the dependence on the coordinate space variables x and y 
and T'^ (T") represent the (anti-)time-ordering operators. In the definition of the fac- 
tor rj = +1 for bosons and rj = -1 for fermions. In the following we consider for simplicity 
a theory for scalar bosons. The modifications for relativistic fermions or nonrelativistic 
particles will be presented in connection with the final equations. The full Green functions 
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are determined via the Dyson-Schwinger equations for path-ordered quantities, here given 
in 2 X 2 matrix representation 



q> qa I I q> qa I ~'~ I q> qa I \ _V> V" / I q> q<^ 

'-'xy '-'xy / \ '-'0,xy '-'0,xy / \ '-'0,xzi '-'0,xzi / \ ^2122 ^2122 / \ 221/ '-'zzy 

The self-energies S' are also defined according to their time structure while the symbol 
"0" implies an integration over the intermediate spacetime coordinates from — oo to oo. 
Linear combinations of diagonal and off-diagonal matrix elements give the retarded and 
advanced Green functions 5'^'^*/"'^^ and self-energies J]^^^/"-''-'" 



qret qc q< q> qa qadv qc q> q< qa 

xy xy xy xy xy ' xy xy xy xy xy ' 

y^ret ^adv v^c 

xy xy xy xy xy ' xy xy xy xy xy " \ ) 

Resorting equations (^) one obtains Dyson-Schwinger equations for the retarded/advanced 
Green functions (where only the respective self-energies are involved) 

Q— 1 Qret,adv r i \^ret,adv ^ Qret,adv f a\ 

'-'Ox '-'xy — ^xy "T l-'^z ^ '~>zy ) 

and the wellknown Kadanoff-Baym equation for the Wightman function S^, 

^Ox ^xy = ^xz ® ^^y + ^xz © ^zy^ ■ (5) 

In these equations 5*1^^ denotes the (negative) Klein-Gordon differential operator which 
for bosonic field quanta of (bare) mass Mq is given by Sq^ = —{d^d^ + Mq). The Klein- 
Gordon equation is solved by the free propagators 5*0 as 

q-l I ^0,xy '^0,xy \ _ X I ^ ^ \ q-l qret,adv _ JT /c\ 

^Ox C> qa j — ^xy \ n. 1 ' ^Ox ^0,xj/ — "xy 

\ '^0,xy '-'0,xy ) \ ^ / 

with the four-dimensional 5-distribution 5xy = 5^^\x — y). 

In the following one changes to the Wigner representation via Fourier transformation of 
the rapidly oscillating relative coordinate (x — y) and formulates the theory in terms of 



the coordinates X = {x + y)/2 and the momentum P, 



FxP = / d\x-y) e'^^^^'-y'^ F^y. (7) 

Since convolution integrals convert under Wigner transformations as 

jd\x-y) e'^^^-'-y') Fi,,,0F2,,, = e^'^ F,,px F2,px, (8) 

one has to deal with an infinite series in the differential operator O which is a four- 
dimensional generalization of the Poisson- bracket, 

0{Fi}{F,} .- ^[g^^Q^, - 9^9^J ■ (9) 

As a standard approximation of kinetic theory only contributions up to first order in the 
gradients are considered. This is justified if the gradients in the mean spatial coordinate 
X are small. 

Applying this approximation scheme (Wigner transformation and neglecting all gradient 
terms of order n > 2) to the Dyson- Schwinger equations of the retarded and advanced 
Green functions one ends up with 

(p2 _ ^2 _ jieJ:^^) ReS^P = 1 - \rxpAxp , 

(p2 _ ^2 _ R^j^r^t^^ AxP = TxpReS^'x% , (10) 

where we have separated the retarded and advanced Green functions as well as the self- 
energies into real and imaginary contributions 

STp''"'' = ReS^P t\Axp, TT^'i^'' = i?eE-*, T ^^xp . (11) 

The imaginary part of the retarded propagator is given (up to a factor) by the normalized 
spectral function 



Axp = i 



oret oadv 



r HP^ 

^-2ImSTp, J^AxP^l, 



5 



while the imaginary part of the self-energy is half the width Txp- From the algebraic 
equations (|iy) we obtain a direct relation of the real and the imaginary part of the 
propagator (provided Txp > 0): 

P2 _ /\/f2 _ Dpyret 

ReSTp = ^'^^^ AxP . (12) 

J- XP 

Finally the solution for the spectral function shows a Lorentzian shape with space-time 
and four-momentum dependent width Txp- This result is vahd for bosons to the first 
order in the gradient expansion, 

(P2 _ Mi - ReH^^'pf + T\p/A ■ ^ ^ 

For the real part of the Green function we get 

r^^^ret ^ - - ReT^^P .... 

(P2 _ Mi - ReT.'^i'pY + V\p/A ' ^ ' 



2.2 Transport equations 

In this Subsection we derive the transport equations that will be used to describe the prop- 
agation of particles with a space-time dependent finite life time. For this aim we start 
with the Kadanoff-Baym equation (|^) which yields in the same approximation scheme 
(i.e. a first order gradient expansion of the Wigner transformed equation) by separating 
real and imaginary parts to a generalized transport equation, 

0{P' - Ml - ReY.%}{S<p} - 0{T.<p}{ReSTp] 

~ 2 [ ^^-P ~ ^xp ^xp ] 5 (15) 

and to a generalized mass-shell equation, 

[P' - Ml - PeS-^]5<p - T.<pReSTp 

= lo{^lp}{Axp} - lo{rxp}{S^p}. (16) 

In the transport equation (|T3|) one recognizes on the l.h.s. the drift term P^(9^«, generated 
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by the contribution 0{P'^ — Mq}{»}, as well as the Vlasov term determined by the real 
part of the retarded self-energy. On the other hand the r.h.s. represents the collision 



term with its 'gain and loss' structure. To evaluate the 0{S^}{i?e5'^'^*}-term in (|T^), 
which does not contribute in the quasiparticle limit, it is useful to introduce distribution 
functions for the Green functions and self-energies as 

1 S<p = Nxp Axp , iS>p = {I + Nxp ) Axp , 

2 E<p = iV|p r^P , zS>p = (1 + iV|p)rxp. (17) 

Following the argumentation of Botermans and Malfliet the distribution functions 
and A^^ in (|17D should be equal within the second term of the l.h.s. of (|15D within a 
consistent first order gradient expansion. In order to demonstrate their argument we write 

S<p = -iTxpN^p = -iTxpNxP + Cxp. (18) 



The 'correction' term Cxp is proportional to the collision term (r.h.s.) of the generalized 
transport equation (|15|), 



Cxp = -tTxpiN^p - Nxp) = z ( S< p 5>p - S>p 5<p ) A^^p , (19) 

which itself is of first order in the gradients. Thus, whenever a distribution function 
appears within a Poisson bracket the difference term (A^^ — A^) becomes of second order 
in the gradients and has to be omitted for consistency. As a consequence A^^ can be 
replaced by A^ and the self-energy by 5^ • T/A in the term 0{S^}{i?eS''''^*}. The 
general transport equation (|15|) then can be written as 



A 



XP J- XP 



0{P'-M^-ReEfp}{S<p} - -L<y{Txp}{{P'-M',^ReETp)S<p} 

-L XP J 

= i [ S^p Sxp — ^xp ^xp ] • (20) 



In order to explore the physical content of ( [20| ) we simplify the problem by assuming the 
self-energy only to depend on space-time coordinates, i.e. ReT,xp = -ReS^*,rxp = Tx- 
In view of transport applications it is now advantageous to introduce a 'virtual mass 
parameter' that is determined by the four-momentum of the particles and incorporates 
also energy-shifts generated by the real part of the self-energy as 

= - ReJ:f . (21) 
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Taking as an independent variable this fixes the energy (for given P and M^) to 

= P^ + + ReTrs\ (22) 
In this limit the transport equation simplifies as 

0{P'-Ml~ReTri'p}{S<p] - ^0{Txp}{{P'-Ml~Rei:^i'p)S<p} 

J- XP 



-^(^^Tx) [d''p[{P'-Ml-Re^f)S<p]) . (23) 



To solve equation (pOD we use a testparticle ansatz for the real quantity p j^.p , which 
represents the probability to find a particle in a given phase-space cell, 

N 

FxPM^ = ^^XPM^ - E - Ut)) 5^'\P - m) 6{M' - Mm) ■ (24) 

i=l 

In ( |2^ ) we have formulated the testparticle ansatz in terms of the 'mass parameter' 
instead of the energy Pq to gain the evolution in M^. F-^pj^^ then is a solution of (|23|) if 
the testparticle coordinates obey the following equations of motion 



dt P 



(25) 



Oi 



dMf Mf - Mj d Mf-Mj P 

For Tx = we obtain directly the familiar equations of motion in the quasiparticle 
approximation assuming quasiparticles states with an effective mass squared Mq + PeS^* 
and a spectral function proportional to a (5-function. This limit has been traditionally 
employed in transport theories H, ^ allthough its applicability should be restricted 
to low collision rates. Starting from the general equation (pOD, which does not invoke 
any on-shell approximation, the quasiparticle limit can also be obtained for small width 
r by integrating out the energy Pq. The second term on the l.h.s. of (113]) then vanishes 



while all other terms are given in terms of corresponding quasiparticle quantities (Green 
functions, self energies, etc.). We note that the quasiparticle approximation is known 



8 



to be a self-consistent energy conserving approximation. Furthermore, the conventional 
quasiparticle approximation also holds for particles of finite life time if dt T x = V Tx = 0, 
i.e. in the low density regime with almost vanishing collision rate Tcoii (see below). 
Eqs. - ( P?!) fulfill energy conservation with respect to the particle propagation if 

^.v.r, -|r, = o, (28) 

i.e. dM'^/dt = —2P-dP/dt (with respect to the additional terms involving Tx)- Eq. ( pSf ) 
is identically fulfilled for dTx/dt = according to (P7|), however, violated for dVx/dt ^ 0. 
In view of a semiclassical treatment we suggest to fix the energy Poi at every time t by 
adding a term ~ dTx/dt in eq. (|26|), i.e. 



dt Po, 



(29) 



^ - r^. rf^^^- ^^^^ 

Dynamical effects of the additional term in (BUp will be discussed in the first applications 



of our model in Section 3. 

In case of relativistic fermions with (bare) mass uiq we have 

^ret T ^2 I T/'2 i oTD oo T/ i o,™ rr i rrS 



PeE^^^* = -Vq' + V + 2PoVo - 2P -V + 2mo [/cj + f/| , (32) 



where = (Vb, V) and f/^ denote the real part of the vector and scalar self-energy |^ 



respectively. Here we have neglected pseudo-scalar, pseudo-vector and tensor contribu- 
tions in the Gordon decomposition of the self-energy, which either vanish identically or 
remain small in the nuclear physics context. 
A similar decomposition holds for Tx, i.e. 

Tx = -W^ + + 2PoWo - 2P -W + 2mo Ws + (33) 

for the vector part (W^) and scalar part (Ws) of the imaginary self-energy. 
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3 Model studies 

For our present purpose to demonstrate the physical imphcations of eqs. ( P^ - ( PT| ) we 
consider the propagation of particles in a time-dependent complex potential of Woods- 
Saxon form, i.e. 



- - { l + exp{(|^°-i?)/ao} " ' [l + e^vm-R)/ao] + t)} ^^^^ 

where we have used R = 5 fm, oq = 0.6 fm throughout the model studies. Eqs. (^) - 
(|3TD allow to represent the distribution function in terms of the testparticle distribution 
( p^ where ri(t),Pi(t) and Mf(t) are the corresponding solutions of eqs. ( pQ]) - (pT]). We 
initialize all testparticles i with a fixed energy Pq at some distance (|r'(t = 0)| ~ —15 fm) 
on the 2;-axis with a three-momentum vector in positive 2;-direction. The mass parameters 
Mi{t = 0) are selected according to the Breit-Wigner-distribution 

F(M) = —- (35) 

^ ^ 271 {M - Mo)2 + r^/4' ^ ' 

where Ty denotes the vacuum width which might be arbitrarily small but finite (see 
below). The particles are then propagated in time according to eqs. (^) - ( |3TD and all 
one-body quantities can be evaluated from (^). 

In Fig. 1 (upper part) the results for Pio{z{t)), Mi{z{t)) and Piz{z{t)) are diplayed as a 
function of z{t) instead of the time t. We show the evolution of 21 testparticles with mass 
parameters that are initially seperated by AM = 0.05 ■ Ty in the case of a nonvanishing 
imaginary part of the potential (Wo{t) = Wq = 70 MeV, Fy = 0.8 MeV) but vanishing 
real part of the potential (Vq = MeV) (see Fig. 1 (lower part)). One recognizes that 
the differences between the mass parameters increase when reaching the potential region, 
which corresponds directly to a broadening of the spectral function. The same spreading 
behavior is observed for the three-momentum of the testparticles, such that the energy Pq 
is conserved throughout the whole calculation (upper line). When leaving the potential 
region the splitting decreases and the correct asymptotic solution is restored. 

At next we present in Fig. 2 (upper part) a calculation where we additionally allow 
for a nonvanishing real part of the potential (i.e. Vq = —20 MeV, Fig. 2 (lower part)). 
While the spreading of the mass parameter is not affected by this change, we find a shift 
of the testparticle momenta where the real part of the potential deviates from zero since 
here the particles are accelerated. 

While we up to now have only considered constant potentials in time, we now introduce 
an exphcit time dependence corresponding to Wo(t) = 100 MeV (1 — 0.02 c/fm-t). As a 
result the spatial reflection symmetry vanishes for Piz{z(t)) and Mi{z{t)) (cf. Fig. 3). 
For the given time dependence the mass splitting is smaller for given z compared to a 
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time-independent potential. As in the former cases here also the correct free solution is 
obtained for large z while the energy is strictly conserved, too. 

In order to discuss the effect of the additional term ~ dVx / dt introduced in eq. ( pOD 
we solve eqs. (^) - (^) that directly stem from the testparticle ansatz for the model 
case presented in Fig. 3. The corresponding results in Fig. 4 show that the time evolu- 
tion in is hardly effected, however, the individual trajectories spread out in energy Pq 
which can be traced back to an asymptotic spread in momentum P^. When integrating 
over the particle spectral function it can be shown that the total energy of the particle - 
which corresponds to a quantum mechanical state - is conserved again due to the sym- 
metry of the spectral function in (M^ — Mq) for a local width F^. Thus the original 
equations of motion ( p5D - ( pT]) conserve energy with respect to the quantum mechani- 
cal state. However, within semiclassical transport simulations the spectral function will 
not be populated symmetrically around Mq due to energy constraints, e.g. in reactions 
ppp^ with invariant energy ^/s < 2Mp + Mq. Thus only the low mass fraction of the 
Breit-Wigner distribution will be populated for which the energy is no longer conserved 
throughout the propagation which might produce artefacts in further inelastic reaction 
channels. Such artefacts do not occur in a fully quantum mechanical theory due to the 
phase coherence ~ iPot which guarantees energy conservation for t ^ oo. Since semi- 
classical transport simulations do not involve such phase coherence and thus violate the 
uncertainty relation with respect to energy and time, we have restored energy conserva- 
tion locally in eqs. (^) - (^) for each testparticle representing a tiny slice in dM"^ of the 
spectral function. 

In the next example of this model study we show in Fig. 5 (upper part) the case of 
a broad vacuum spectral function entering a (time- independent) nonrelativistic potential 
with Vo = -20 MeV and Woit) = Wq = 100 MeV. The vacuum width is chosen as 
Ty = 160 MeV, while 11 testparticle trajectories are shown with an initial separation of 
the masses AM = 0.05 ■ Ty- One observes that the spectral function is further broadened 
in the complex potential zone and reaches its inital dispersion in mass again after passing 
the diffractive and absorptive area. 

The question remains if the testparticle distribution ( p4D reproduces the local splitting 
in mass as expected due to quantum mechanics, i.e. in our case a Breit-Wigner distribution 
(|35D with a local width Fx = 2 Wo{z) + Ty- This is demonstrated in Fig. 6 where we 
show the spectral function as a function of mass M from the testparticle evolution at fixed 
coordinate z in comparison to the quantum Breit-Wigner distribution with local width 
Tx (full lines) for a pure imaginary potential with parameters Wq = 50 MeV and vacuum 
width Fy = 2 MeV. The differences from the exact results in Fig. 6 are practically not 
visible for all values of z from - 8 fm to 8 fm. The width of the distribution increases 
from 1 MeV in the vacuum {z = ± 8 fm) to 102 MeV (= 2 Wq+ Fy) in the center of the 
absorptive potential {z = 0). Thus our off-shell quasiparticle propagation is fully in line 
with the quantum mechanical result at least for quasi-stationary quantum states. 

To summarize our model results for the simple complex potential of Woods-Saxon- 
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type, we find that energy conservation is guaranteed during the propagation as well as 
that the correct asymptotic solutions for the spectral functions are restored. Furthermore, 
in the potential region we observe a broadening of the width of the spectral function due to 
the space-time dependent imaginary part of the potential in line with quantum mechanics. 



4 Application to nucleus-nucleus collisions 

Apart from the model studies performed in the previous Section it is of interest, if the 'off- 
mass-shell approach' proposed here leads to observable consequences in actual experiments 
such as for nucleus- nucleus collisions. This implies to specify the collision term in eq. (p!5|). 
A corresponding expression can be formulated in full analogy to Ref. by giving explicit 
approximations for and and using detailed balance as 

hou{X, P, = Tr2TrsTuA{X, P, M^)A{X, P^, M|) A(X, P^, Mi)A{X, P4, M^) (36) 
{|T((P, M^) + (P2, ^ (P3, Mi) + (P4, Ml))\Xs S\P + P2-P3- P4) 

[^XP3M^^XPiMlfxP,kpfxP2,Ml ~ ^ XPM'^^ XP2M^f XP-iMs^ XPiMl^^ 

with 

f^P = l + vN^P (37) 

and 1] = ±1 for bosons/fermions, respectively. The indices A, S stand for the antisymmet- 
ric/symmetric matrix element of the scattering amplitude T in case of fermions/bosons. 
We note that in (|36|) we have neglected gradient terms which lead to nonlocal collision 
terms ^ Their effect can additionally be taken into account following 



the suggestions of Ref. , but are discarded here for simplicity and transparency. 
In eq. (^) the trace over particles 2,3,4 reads explicitly for fermions 

^ 1 /• ,0 dM^ 



where ct2, T2 denote the spin and isospin of particle 2. In case of bosons we have 

Tr^ = V / d'P2 (39) 

since here the spectral function is normalized as 
r/p2 

-^As{X,P) = l (40) 
whereas for fermions we have 
lfMX,P)^^^. (41) 



12 



It is easy to show that the colhsion term (^) leads to the proper Bose or Fermion equi- 
hbrium distributions for t ^ oo. 

Neglecting the 'gain-term' in eq. ( ^61) one recognizes that the coUisional width of the 
particle in the rest frame is given by 

T,oii{X, P, M^) = Tr^Tr^Tu {|T((P, M^) + (P2, Mj) ^ (P3, Ml) + (P4, Ml))\\s (42) 

A(X, P2, M|)A(X, P3, Ml)A{X, P,, Ml)5\P + P2-P3- P4) N^P.MifxpJxpJ^ 

where as in eq. (^) local on-shell scattering processes are assumed. We note that 
the extension of eq. (^6]) to inelastic scattering processes (e.g. A^A^ — A^A) or {ttN 
A ect.) is straightforward when exchanging the elastic transition amplitude T by the 
corresponding inelastic one and taking care of Pauli-blocking or Bose-enhancement for 
the particles in the final state. We mention that for bosons we will neglect a Bose- 
enhancement since their actual phase-space density is small for the systems of interest. 
For particles of infinite life time in vacuum - such as protons - the collisional width 
has to be identified with half the imaginary part of the self-energy that determines 
the spectral function (|35|). Thus the transport approach determines the particle spectral 



function dynamically via (52) for all hadrons if the in-medium transition amplitudes T 
are known. Since in binary collisions due to energy and momentum conservation - once 
the final masses are fixed - only the final scattering angle fl = (cos 6*, 0) is undetermined 
we can replace the amplitude squared in (^2]) as 

where q is the momentum transfered in the collision at invariant energy ^/s and /i is the 
reduced mass of the scattering particles. The differential cross section da/dQ or T{q) in 
principle should be evaluated in the Brueckner approach, however, in practice effective 
parametrizations are employed (see below). 



4.1 Numerical realisation 



The following dynamical calculations are based on the conventional HSD transport ap- 
proach [13, ^], where for energies up to 100 A MeV (GANIL energies) essentially the 
nucleon degrees of freedom are important, since inelastic processes A^A^ NA 
7tN, 7cN A are suppressed. We only briefly mention that the formation cross sec- 



tion in the reaction A^A^ 



A^A is evaluated with the total A-width F^* 



^ A 



+ F 



coll 



41| and that in the decay channel A — > vrA^ the final nucleon state is selected by Monte 



Carlo using the local spectral distribution 

Whereas the real part of the nucleon self-energy is determined as in Ref. 
an explicit momentum dependence of the scalar and vector self-energies for nucleons in 



131 and includes 
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order to qualify also for relativistic reactions, we have to describe in more detail the im- 
plementation of the off-shell dynamics induced by Txp = Tcoii{X, P, M"^). According to 
(|42D the coUisional width is explicitly momentum (and energy) dependent, which intro- 
duces much larger numerical efforts as for momentum-dependent real potentials. In view 
of the limited energy range addressed here and for the purpose of an exploratory study it 
is sufficient to consider a space-time dependent collisional width Fx which is obtained by 



averaging (^) at each time-step (At = 0.5/m/c) in each cell in coordinate space (of size 

with the local density 

1 r rlM"^ 

We note, that in order to achieve a numerically 'flat' function Tx in space and time, 
one has to consider averages over a large set of ensembles typically in the order of 10^ 
testparticles per nucleon. By storing Tx on a 4-dimensional grid the space-time derivatives 



of Txi that enter the equations of motion (^9]) - (0), can be evaluated in first or second 



order. We mention that the Gaussian smearing algorithm described in Ref. leads to 
sufficiently stable results (see below). 

The collisions of nucleons are described by the closest distance criterion of Kodama et al. 
l5l| in the individual A^A^ c.m.s., i.e. 



|Xi-X2|<V^(^i-A)/7r (46) 
using the Cugnon parametrizations |E2| for the in-medium A^A^ cross section da/dQ{y/s) 



and identifying (in the A^A^ c.m.s.) 

s - Am% = s- 4M2 = 4P^ (47) 

where is the nucleon vacuum mass, M the actual off-shell mass and y/s the invariant 
energy of a nucleon-nucleon collision in the vacuum with c.m.s. momentum P. 



According to eq. (|36D the nucleons can change their virtual mass M in the scattering 
process 1 + 2-^3 + 4, while keeping the energy and momentum balance. This process is 
technically handeled by selecting the final nucleon masses by Monte Carlo according to 
the local Breit-Wigner distribution. However, our Monte Carlo simulations showed that 
this change of virtuality for elastic collisions has a minor effect on the observables to be 
discussed below. 

Apart from the description of particle propagation and rescattering the results of the 
transport approach also depend on the initial conditions, Xj(0), Pi(0), Mp(0). In view of 
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nucleus-nucleus collisions, i.e. two nuclei impinging towards each other with a labora- 
tory momentum per particle Piab/A, the nuclei can be considered as in their respective 
groundstate, which in the semiclassical limit is given by the local Thomas-Fermi distribu- 
tion 0. Additionally the virtual mass has been determined by Monte-Carlo according 
to the Breit-Wigner distribution (|35|) assuming an in-medium width Fq = 1 MeV. For the 
vacuum width of the nucleons we have used Fy = 1 MeV which implies that nucleons 
propagating to the continuum in the final state of the reaction achieve their vacuum mass 
on the 0.1 % level. We note in passing that for the initialization we additionally have 
required Pqi < "^tv (in the restframe of the nucleus) which due to energy conservation 
implies that particles cannot escape 'numerically' from the nucleus in the groundstate. 



4.2 Nucleus-nucleus collisions at GANIL energies 

Our first applications we devote to nuclear reactions at GANIL (92 - 95 A MeV) since here 
the more recent measurements have lead to confiicting results between different transport 



approaches |£3[. We start with the reaction Ar + Ta at 92 A MeV. 
In view of Section 3 we present for some randomly chosen testparticles i their off-mass- 
shell behaviour Mf{t)-M^ function of time in a central collision (6=1 fm) in Fig. 7. 
It is seen that during the maximum overlap of the nuclei at t ~ 30 fm/c the off-shellness 
reaches up to 0.2 GeV^, however, in analogy to the model studies in Section 3 the nucleons 
become practically on-shell for t > 90 fm/c. The finite width at the end of the calculation 
presented here is due to the fact that the coUisional width Tcou is still different from zero. 
The fluctuations in Mf{t) — Mq in time give some idea about the numerical accuracy of 
the calculation for the space-time derivative of Tx'-, the functions become smoother when 
increasing the number of testparticles/nucleon furtheron (> 1000). 

Without explicit representation we note that the proton rapidity spectra dN/dy do not 
change within the numerical accuracy when comparing the on-shell propagation limit 
with the results from the off-shell transport approach for Ar + Ta at 92 A MeV. There is, 
however, a small enhancement in the proton transverse momentum spectra l/prdN/dpT 
for the off-shell propagation of nucleons as can be seen in Fig. 8, when the proton pt 
spectra are compared at an impact parameter 6=1 fm. 

The question remains, if such an enhancement might be seen experimentally or if the 
off-shell approach overpredicts the high momentum tail of the spectra. For this purpose 
we compare to the proton spectra from Ref. taken at 6iab = 75° for kinetic energies 



above 175 MeV (Fig. 9). For this comparison we have integrated the proton spectra over 
all impact parameters with stepsize Ab = 1 fm in the angular range 70° < 6iab < 80°. The 
calculated spectra (open squares in Fig. 9) only extend up to 200 MeV due to statistics, 
however, match with the experimental data within the errorbars. A similar comparison 
has been performed by Germain et al. |5^ where the traditional BUU and QMD calcu- 
lations seem to be compatible with the experimental data within the statistics achieved 
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whereas the Boltzmann Langevin (BL) calculations incorporating fluctuations in momen- 
tum space overestimate the proton spectra by more than three orders of magnitude 
[p3| . This rules out the latter BL approach but does not imply that our off-shell transport 
approach will properly describe the experimental spectra up to Ekin = 350 MeV. In view 
of Fig. 9 the calculational statistics would have to be increased by more than 3 orders of 
magnitude which is unlikely to achieve for our present off-shell approach within a reason- 
able time due to the high amount of processor capacity needed. Note that due to nonlocal 
effects in the collision term the proton spectra might be slightly hardened additionally 



m 



We thus continue with more qualitative investigations, that allow to extract the physics 
more clearly. In Fig. 10 we display the number of baryon-baryon collisions dN^^ / d^/s 
as a function of the invariant energy ^/s for Ar + Ta at 92 A MeV integrated over all 
impact parameters. The dashed line shows the result for the on-shell transport approach 
(starting at 2mi\f) whereas the solid line corresponds to the off-shell result, which extends 
down to ^/s ^1.5 GeV. Note that elastic collisions of off-shell nucleons can occur due to 
their dynamical virtuality in mass. The dashed line is practically identical to the BUU 
calculations from Ref. and is limited to collisions far below the kaon production 
threshold of y/so ~ 2.54 GeV. Thus the kaon production yield of (2.9±1.6) ■ 10~^b claimed 



in Ref. |^ for this system cannot be described in the on-shell limit, however, also not in 
our off-shell approach which shows only a small enhancement in the high ^/s regime. 
The latter ^/s distribution can approximately be tested experimentally by hard photon 
spectra, a question that has been explored by the TAPS collaboration for Ar + Au at 95 



A MeV |5J] . In order to test our transport approach we have performed calculations 
for this system, too, using the parametrizations (4.13) of Ref. [^j for the elementary dif- 
ferential photon cross section in proton-neutron {pn) collisions. Note that the elementary 
photon bremsstrahlung in pn collisions is at best known within a factor of 2 (cf. the 
discussion in Ref. 0). Fig. 11 displays the results of our bremsstrahlung calculations 
in comparison to the data from Ref. The dashed line corresponds to the conven- 



tional on-shell calculation and is practically identical to the BUU analysis performed by 
Holzmann et al. |3 , but underestimates the high energy photon yield dramatically. This 
situation does not improve very much when including the off-shell propagation of nucleons 
for the initial channel (dash-dotted line), however, still requiring that the nucleons in the 
final state are on-shell, too. Denoting off-shell nucleons by an extra * this corresponds to 
the individual reactions p* + n* —>■ p + n + 'j, whereas the dashed line is obtained from 
the channel p + n—yp + n + 'j. 

On the other hand, the energetic photons are produced very early in the collision phase 
where the virtual mass distribution of nucleons - determined by Tx (^21) - becomes very 
broad. Thus including this virtuality in mass also in the final state, where the masses 
are selected by Monte-Carlo according to (|35D with a local width Tx, we selfconsistently 
can sum the individual channels p* + n* ^ p'* + n'* + 7. The result of such calculations 
is shown in Fig. 11 by the solid line which comes quite close to the experimental data 
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47| . We note that in the latter calculations we have averaged the photon yield over 



10 MeV bins to reduce the statistical fluctuations emerging from the Monte-Carlo final 
state selection. Whereas in Ref. 0] the high energy photon yield has been tentatively 
attributed to very high momentum components in the initial phase-space distribution - 
which semiclassically are not bound - our present results indicate that this yield might 
be almost entirely explained (without introducing any additional assumptions) by the 
off-shell transport approach. It is presently unclear, if the missing high energy photon 
yield should be attributed to three-body reaction channels to the contribution of 

the A 'yN channel |57] or to the secondary vrA^ 'yN channel [5S, 59 . 
The question now arises if the kaon yield from Ref. for Ar + Ta at 92 A MeV might 
also be due to off-shell hadronic states in the final channel. We thus have performed cal- 
culations for this system again within the following assumptions: = Tj^ and = 
since the KN cross section is rather small 



eg 



The parametrisations for the elementary 
reactions A^A^ —>■ NAK and vrA^ —>■ KA have been adopted from Refs. 61] where 
K, K production has been systematically investigated for nucleus-nucleus reactions in the 
SIS energy regime. In view of Fig. 11 the small enhancement of energetic collisions for 
off-shell nucleons in the entrance channel does not have a large effect on the kaon produc- 
tion cross section as in case of energetic photons (cf. Fig. 11), but the reduction of the 
threshold due to the virtuality in mass of the final nucleon and A hyperon should have 
(cf. Fig. 11). Unfortunately, within the statistics achieved in our calculations we did not 
find any KA production event such that we presently cannot provide a solid answer. 
In order to obtain some upper estimate we parametrize the coUisional distribution dN^"-/ dy/s 
(cf. Fig. 10) for pn collisions by an exponential tail 



d^s Eq 



(48) 



with some slope parameter Eq and fix the normalization constant as well as Eq by the 7 
spectra from the TAPS collaboration. The result of this model study is shown in Fig. 12 
where the computed photon spectra for £^0 = 33 MeV (dashed line) and Eq = 35 MeV 



(solid line) are shown in comparison to the data We note that within this model the 
photon spectra remain roughly exponential up to 700 MeV photon energy, which would 
have to be proven by experiment explicitly. 

Now replacing the elementary differential photon cross section dcTpn^pn-yi^/s) / dVt by the 
elementary kaon production cross section daiyN^KAN{\/s) |60| and correcting for isospin 
we obtain an inclusive kaon production cross section ~ 1.6 ■ lO^^^b {Eq = 33 MeV) 
and ~ 6 ■ 10~^^6 {Eq = 35 MeV) for Ar + Ta at 92 A MeV, which is about two orders 
of magnitude smaller than the cross section from Ref. 
has to be added here: In view of the analysis in Refs. 



48 



ra, 62 



Moreover, a note of caution 
i, 1611, |63, |6l the kaon 



potential in the nuclear medium is most likely repulsive and the A potential only 2/3 of 
the nucleon potential. This shift in production threshold to higher y/s due to the real 
part of the hadron self-energies makes kaon production even more unlikely. We thus do 
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not expect to describe the kaon cross section from Ref. in our off-shell approach even 
when increasing the statistics by some orders of magnitude. 



5 Summary 

In this work we have derived and developed a semi classical transport approach that in 
first order in the gradient expansion describes the virtual propagation of particles in the 
invariant mass squared besides the conventional propagation in the mean-field po- 
tential given by the real part of the self-energy. The derivation has been based on the 
familiar Kadanoff-Baym equations by exploiting the relations between the different 



Green functions and the relations between their real and imaginary parts. Whereas in 
conventional transport approaches the imaginary part of the self-energy is reformulated 
in terms of a collision integral and simulated by on-shell binary collisions, we addition- 
ally account for the off-shell propagation of particles due to the imaginary part of the 
self-energy in eqs. ( ^Of ) - (|3lD . We note that in our formulation the single-particle energy 



Pq is fixed by eq. ( p2D and that the propagation is determined by dP^/dt = in the 
coUisionless limit (if ReJ^x^ has no explicit time dependence). On the other hand, the 
local collision rate Tx is determined by the collision integrals themselves and can be used 
in transport approaches without introducing any new assumptions or parameters. In our 
present approach we have restricted to momentum-independent imaginary self-energies; 
an extension to the general case appears straight- forward according to Section 2, however, 
is numerically much more involved. 

As a first application we have studied the dynamical evolution of particles in a fixed 
complex potential - having some similarities to hadron-nucleus collisions without explicit 
collisions - and demonstrated the off-mass- shell propagation in a transparent way for a 
variety of model cases. As also shown numerically, the energy conservation strictly holds 



for the set of eqs. (^) - (0). Furthermore, a distribution of off-shell particles regains 
its vacuum spectral function when moving out of the complex scattering centre; particles 
with vanishing (or very small) width become asymptotically on-shell again as required by 
the quantum mechanical boundary conditions while in the absorptive medium the spec- 
tral function reproduces the correct width in line with quantum mechanics (cf. Fig. 6). 
We have, furthermore, presented the first dynamical calculations of the novel transport 
theory for nucleus-nucleus collisions at GANIL energies where we can test its results in 
comparison to experimental data. We find that the off-shell propagation of nucleons 
practically does not change the rapidity distributions dN/dy and only has a moderate 
effect on the high transverse momentum spectra of protons. The latter we found to be 
fully compatible with the data from Ref. |^ for Ar + Ta at 92 A MeV contrary to the 



Boltzmann-Langevin (BL) calculations in Ref. ||53[. The distribution of nucleon-nucleon 



collisions in the invariant energy ^/s is found to be also slightly enhanced for high in- 
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variant energies (as well as below the 2 nucleon threshold), which has some effect on the 
production of high energy 7-rays. Here we controlled our calculations by the photon data 
from the TAPS collaboration for Ar + Au at 95 A MeV [0 that could be reasonably 



described in the off-shell limit when including especially the nucleon spectral functions 
in the exit channel p* + n* + 7. We have argued that this might be a first experimental 
indication for the off-shell propagation of nucleons that has to occur due to quantum 
mechanics. Our attempt to calculate the kaon cross section for Ar + Ta at 92 A MeV |^ 



within the same line failed due to the limited statistics. However, when extrapolating the 
coUisional -y/i distribution by an exponentail tail ( ^HD and fixing the slope parameter by 
the photon data from the TAPS collaboration our upper limit for this cross section 
is roughly two orders of magnitude below the experimental value, which implies that the 
data point |^ remains ununderstood theoretically furtheron. 

The actual applications of the present off-shell transport approach are not limited to nu- 
clear physics problems, but should be of relevance for any system including particles of 
finite life time and/or high collision rates. The more practical point is now to set up 
new numerical recipies to increase statistics and also to include the explicit momentum 
dependence stemming from the imaginary part of the particle self-energy. 



The authors like to acknowledge stimulating discussions with C. Greiner and S. Leupold 
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ment. 
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Figure 1: upper part: Pjo, Mi and Piz as a function of z{t) for a purely imaginary potential 
Wo(t) = Vl^o = 70 MeV (lower part). The vacuum width is Fy = 0.8 MeV and the initial 
separation in mass of the testparticles is given by AM = 0.05 • Ty- 
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Figure 2: upper part: Pio, Mi and Piz as a function of z{t) for a complex potential with 
Vo = -20 MeV, Wo{t) = Wo = 70 MeV (lower part). For the vacuum width and the 
initial mass separation we have used the same values as in Fig. 1. 
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Figure 3: upper part: Pjo, and Piz as a function of z{fy for an explicitly time-dependent 
imaginary part of the potential W^it) — 100 MeV (1 — 0.02 c/fm-i). The real part is taken 
as in Fig. 2 (lower part). For the vacuum width and the initial mass seperation we have 
used the same values as in Fig. 1. 
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Figure 4: Pio, Mi and Piz as a function of z(t) for the same time-dependent imaginary- 
part of the potential Wo(t) as in Fig. 3 but discarding the term ~ dVx/dt in eq. (|30D . 
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Figure 5: upper part: Pio, Mi and Piz as a function of z{t) for a broad vacuum spectral 
function in a time-independent potential Vq — —20 MeV, Wo{t) — Wo — 100 MeV (lower 
part). We have chosen a vacuum width Ty = 160 MeV and an initial mass seperation of 
AM — 0.05 • Ty for the testparticle trajectories displayed. 
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Figure 6: The spectral distribution at different coordinates z from the testparticle dis- 
tribution in comparison to the analytical result (solid lines) for Vq = 0, Wq — 50 MeV 
and Ty = 2 MeV. The analytical result is practically identical to the histograms from the 
testparticle distribution and thus hardly visible. 
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Figure 7: The off-shell propagation in mass {Mf{t) — Mq) as a function of time for 16 
randomly chosen testparticles. The system is Ar + Ta at 92 A MeV and impact parameter 
6=1 fm. 
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Figure 8: The transverse momentum spectra of protons for Ar + Ta at 92 A MeV and 
impact parameter 6=1 fm. The dashed histogram is the result from the on-shell propaga- 
tion while the solid histogram is obtained including the off-shell propagation of nucleons. 
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Figure 9: The inclusive kinetic energy spectra of protons for Ar + Ta at 92 A MeV 
including the off-shell propagation in the transport approach (open squares) for 70*^ < 
diab < 80*^ in comparison to the experimental data from |^ (full circles). 
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Figure 10: The number of baryon-baryon (BB) collisions as a function of the invariant 
energy ^/s for Ar + Ta at 92 A MeV integrated over all impact parameter. The solid line 
is obtained from including the off-shell propagation in the transport approach while the 
dashed line stands for the result in the on-shell limit. 
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Figure 11: The inclusive differential photon spectra for Ar + Au at 95 A MeV within 



various limits in comparison to the data from Ref. [47|. The dashed line is obtained in 
the on-shell propagation limit including on-shell nucleons in the final state, too. The dash- 
dotted line results from the off-shell propagation, however, including on-shell nucleons in 
the final production channel. The solid line results from the off-shell propagation of 
nucleons including also off-shell nucleons in the final channel. 
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Figure 12: The inclusive differential photon spectra for Ar + Au at 95 A MeV using the 
exponential y/s distribution (^8]) for Eq = 33 MeV (dashed line) and Eq = 35 MeV (solid 
line) in comparison to the data from Ref. [^. 
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